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This  Instruction  Report  was  prepared  at  the  Coastal  Engineering  Research 
Center  (CERC)  of  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  under 
Civil  Works  Research  Work  Unit  31592,  "Wave  Estimation  for  Design,"  Coastal 
Flooding  and  Storm  Protection  Program,  Coastal  Engineering  Area  of  Civil  Works 
Research  and  Development,  Office,  Chief  of  Engineers  (OCE).  The  Technical 
Monitor  from  OCE  was  Mr.  John  H.  Lockhart,  Jr. 

The  program  described  in  this  report  is  a  numerical  model  that  simulates 
shallow-water  wave  growth,  propagation,  and  decay  in  a  directional  spectrum 
over  an  arbitrary  bathymetry.  The  model  was  first  presented  to  US  Army  Corps 
of  Engineers  District  and  Division  personnel  during  a  workshop  held  at  WES  in 
January  1985.  Subsequent  enhancements  to  the  program  have  altered  the  input 
and  output  files  since  the  workshop.  The  model  features  described  in  this  re¬ 
port  are  those  included  in  the  model  at  this  writing.  Future  erhancements  and 
auxiliary  supporting  programs  will  be  documented  as  appendixes  to  this  loose- 
bound  report. 

This  report  was  prepared  by  Dr.  Steven  A.  Hughes,  Research  Hydraulic  En¬ 
gineer,  and  Dr.  Robert  E.  Jensen,  Research  Hydraulic  Engineer,  Coastal  Ocean¬ 
ography  Branch  (COB),  CERC,  under  the  supervision  of  Dr.  Edward  F.  Thompson, 
Chief,  COB,  Mr.  H.  Lee  Butler,  Chief,  Research  Division,  Mr.  Charles  C. 
Calhoun,  Jr.,  Assistant  Chief,  and  Dr.  James  R.  Houston,  Chief,  CERC. 

Dr.  Charles  L.  Vincent  contributed  to  this  report  both  in  the  preparation  and 
in  the  technical  review.  This  report  was  edited  by  Ms.  Shirley  A.  J.  Hanshaw, 
Publications  and  Graphic  Arts  Division,  WES. 

Director  of  WES  was  COL  Allen  F.  Grum,  USA.  Technical  Director  was 
Dr.  Robert  W.  Whalin. 
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A  USER'S  GUIDE  TO  SHAUN:  NUMERICAL  MODEL  FOR  SIMULATION 
OF  SHALLOW-WATER  WAVE  GROWTH,  PROPAGATION,  AND  DECAY 


PART  I:  MODEL  OVERVIEW 
Introduction 

1.  SHALWV  is  a  numerical  model  that  simulates  shallow-water  wave 
growth,  propagation,  and  decay  in  a  directional  spectrum  over  an  arbitrary  ba¬ 
thymetry.  It  is  a  time-dependent  model  requiring  a  primary  input  of  the  wind 
speed  and  wind  direction  at  each  grid  point. 

2.  The  purpose  of  this  Instruction  Report  is  to  provide  an  overview  of 
SHALWV  followed  by  instructions  on  how  to  set  up  the  necessary  input  files  and 
how  to  execute  the  model  on  the  Control  Data  Corporation's  (vendor  for  US  Army 
Corps  of  Engineers)  Cybernet  Computer  System.  The  overview  of  the  model  is 
contained  in  Part  I,  which  gives  the  model  description,  the  theoretical  back¬ 
ground,  the  model  applications,  and  the  model  limitations.  Part  II  provides 
the  information  necessary  to  execute  the  model.  This  includes  model  conven¬ 
tions,  model  input,  and  model  output  with  each  being  illustrated  by  examples. 
Examples  of  the  Job  control  language  (JCL)  specific  to  the  Cybernet  Computer 
System  are  provided  in  Appendix  A. 

Background 

3.  The  numerical  wave  model  SHALWV  was  developed  originally  by 

Dr.  Donald  Resio,  Offshore  and  Coastal  Technologies,  Inc.,  under  contract  to 
the  Coastal  Engineering  Research  Center  (CERC).  The  model  was  delivered  April 
1984,  after  which  work  began  in  moving  the  FORTRAN  source  code  to  the  CYBERNET 
system  and  making  the  changes  necessary  to  compile  and  execute  the  model  on 
the  Cyber  20^  computer.  The  model,  which  was  extensively  enhanced  and  ex¬ 
panded  in  its  capabilities  by  Dr.  C.  L.  Vincent  and  Dr.  R.  E.  Jensen  of  CERC, 
includes  modification  of  the  growth  parameters  to  the  TMA  spectral  shape*and 
the  addition  of  the  ability  to  treat  swell  and  rapidly  changing  winds.  These 
enhancements  have  extended  the  FORTRAN  code  to  about  3,300  lines  in  length. 
Model  Description 

4.  SHALWV  simulates  the  physical  mechanisms  acting  on  a  directional 

*  Steven  A.  Huphec .  1986.  "The  TMA  Shallow-Water  Spectrum:  Decor. rti-  and 

Applications,"  Technical  Report  CERC-fl4-7,  US  Army  Engineer  Waterways  Ex:  - -  - 
mont  Station,  Vicksburp,  Miss. 


spectrum  within  an  arbitrary  depth  field.  The  modeled  spectra  are  represented 
as  fully  two-dimensional  spectra  in  discrete  frequency  and  direction  bands. 
Physical  processes  which  are  modeled  over  an  arbitrary  bathymetry  art-  as 
fol lows : 

a.  Wave  propagation. 

b.  Wave  refraction. 

c.  Wave  shoaling. 

d.  Wave  growth  and  decay . 

These  processes  are  modeled  in  a  time-dependent  mode;  thus,  most  nor.steady- 
state  wave  and  wind  conditions  can  be  handled.  The  model  operates  or,  a  square 
grid  with  directional  spectra  containing  up  to  20  frequency  bands  and  up  to  18 
directional  bands.  The  maximum  time-step  interval  is  determined  from  the  fol¬ 
lowing  Courant  number  stability  criterion  which  is  based  on  grid  size  and 
water  depth: 


At 


>cg(f,) 
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where 


A l  -  grid  spacing 


At  --  time-step 
C  ;  wave  group  speed 

f*  =■  lowest  frequency  band  specified  for  model  during  input 
The  criterion  assures  that  wave  energy  cannot  propagate  more  than  ore  grid 
cell  during  a  t.me-step.  Note  that  C  varies  with  depth. 

O 

5.  Recent  developments  on  depth-limited  spectra  have  been  incorporated 
into  the  model  to  provide  an  upper  limitation  to  wave  growth. 


Computer  r equ  i  reme r. t s 

6.  SHALWV  is  currently  installed  on  the  Cybernet  Computer  System  un¬ 
der  the  family  name  KOE.  Execution  of  the  model  is  through  a  batch  ;oo  that 
transfers  the  necessary  files  to  the  Cyber  205,  which  is  a  large  capacity  "su¬ 
per  computer."  The  Cyber  205  is  used  for  two  reasons:  (a)  the  modei  requires 
a  substantial  amount  of  computer  core  memory  (more  than  is  allowable  on  the 
smaller  interactive  Cybernet  computers',  and  (b)  the  time-dependency  of  the 
model  can  lead  to  leng’  /  computer  runs  which  are  mere  efficient  an  the  much 
faster  Cyber  205.  The  model  conceivably  could  be  installed  on  a  smaller  ma¬ 
chine  with  a  Large  virtual  memory  capacity,  but  this  probably  is  not  a  simple 
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task.  Some  sections  of  the  model  might  have  to  be  modified  to  accommodate  the 
transfer,  and  the  model  execution  would  probably  tie  up  the  entire  cfore  memory 
for  a  long  time.  Therefore,  it  is  recommended  that  the  model  not  be  trans¬ 
ferred  to  a  smaller  machine  except  by  very  experienced  personnel. 

Simulation  costs 

7.  Because  of  the  many  variables  involved  in  model  simulations,  it  is 
possible  to  provide  only  approximate  cost  estimates  for  execution  of  this 
model.  By  way  of  example,  a  12  by  19  grid  of  Lake  Erie  containing  99  water 
grid  points  and  run  for  12  time-steps  (9-hr  simulation)  costs  $18.  The  same 
problem  run  for  29  time-steps  (8  hr)  costs  $23.  A  19  by  15  (92-  by  95-km) 
grid  in  the  Atlantic  off  Duck,  North  Carolina,  ran  for  126  time-steps  (7-hr 
simulation)  at  a  cost  of  approximately  $120. 

Theoretical  Foundation 

8.  This  section  describes  the  theoretical  background  of  SHALWV.  No  at¬ 
tempt  is  made  to  discuss  the  actual  mathematical  formulations  and  the  method 
of  representing  these  formulations  in  the  model. 

9.  SHALWV  employs  an  explicit  finite-difference  scheme  to  solve  the 
following  spectral  balance  equation: 


■?** 


In  this  equation  the  left-hand  term  is  the  total  derivative  of  the  directional 
energy  density,  and  it  contains:  (a)  a  term  representing  the  temporal  change 
in  the  spectral  energy  density,  (b)  terms  representing  the  advection  of  wave 
energy,  and  (c)  refraction  and  shoaling  terms.  The  right-hand  side  of  Equa¬ 
tion  2  consists  of  the  source  terms  representing  physical  mechanisms  that  add 
to  or  subtract  from  net  energy  to  the  wave  field.  These  source  terms  include 
atmospheric  energy  input,  nonlinear  wave-wave  interactions,  bottom  friction, 
bottom  percolation,  wave  dissipation  within  the  wave  field,  bottom  motion, 
bottom  scattering,  and  surf  zone  breaking.  The  numerical  model  SHALWV  con¬ 
tains  the  first  four  source  terms  (atmospheric  input,  wave-wave  interactions, 
bottom  friction,  and  percolation)  and  has  a  finite-depth  spectral  limit  which 
acts  to  constrain  wave  growth  in  lieu  of  a  source  term  for  wave  dissipation  by 
break  Ing . 
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10.  Propagation  of  the  wave  energy  in  each  discrete  frequency-direction 
band  is  done  independently  using  an  upstream  differencing  scheme.  This  is  a 
stepwise  solution  that  estimates  for  each  discrete  band  the  wave  ray  along 
which  the  energy  contained  in  that  band  must  propagate  in  order  to  arrive  at 
the  grid  point  of  interest  by  the  end  of  the  time-step  (Figure  1).  Point  A  on 
Figure  1  represents  the  location  of  the  wave  energy  contained  in  the  discrete 
frequency-direction  band  at  the  beginning  of  the  time-step.  An  estimate  of 
this  energy  is  obtained  by  an  interpolation  method  that  first  projects  the 
wave  ray  farther  back  in  time  until  a  grid  boundary  is  crossed  (point  B  in 
Figure  1).  Finally  this  estimated  energy  is  propagated  along  the  wave  ray  as 
refraction  and  shoaling  effects  are  estimated,  reaching  the  grid  point  (C)  at 
the  end  of  the  time-step.  During  the  propagation,  energy  is  added  to  or  re¬ 
moved  from  each  discrete  energy  band  by  the  source  terms.  At  the  end  of  the 
time-step,  the  directional  spectrum  at  each  grid  point  is  the  sum  of  the  inde¬ 
pendently  propagated  spectral  elements.  Since  the  grid  spacing  can  be  large, 
assumptions  about  smooth  bathymetry  are  inherent. 

Source  terms 

11.  The  SHALWV  model  assumes  that  the  primary  source  mechanisms  for 


Figure  1.  Interpolation  scheme  during  wave 
propagation  over  time-step 


storm  waves  in  shallow  water  are  the  atmospheric  energy  input,  the  nonlinear 
wave-wave  interactions,  and  internal  wave  field  dissipation.  It  is  assumed 
that  bottom  friction  and  percolation  become  important  only  in  certain  in¬ 
stances,  such  as  swell  on  a  gentle  bottom  slope.  Bottom  friction  and  percola¬ 
tion  can  be  activated  by  the  user  if  desired  for  the  swell  components  of  the 
model . 

12.  Figure  2  illustrates  the  relative  effect  that  the  primary  source 
terms  (as  assumed  by  this  model)  have  on  a  wind  sea  spectrum.  The  atmospheric 
source  term  is  formulated  such  that  the  energy  addition  due  to  wind  stress  is 
confined  primarily  at  the  peak  of  the  spectrum  and  to  the  rear  face  of  the 
spectrum  (regions  II  and  III  in  Figure  2)  in  deep  water.  During  fully  devel¬ 
oped  conditions,  the  waves  at  frequencies  less  than  the  peak  travel  at  speeds 
equal  to  or  greater  than  the  wind  speed;  thus  wind  energy  cannot  be  trans¬ 
ferred  directly  to  the  sea  surface.  As  the  water  depth  decreases,  the  lower 
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frequency  waves  become  nond ispersive  and  begin  to  travel  at  a  celerity  pri¬ 
marily  dependent  .toon  the  depth.  SHALWV  assumes  that  as  waves  approach  the 
nond ispersive  limit  in  shallow  water,  the  forward  face  of  the  spectrum  be¬ 
comes  incapable  of  receiving  direct  atmospheric  input,  as  well  as  becoming 
incapable  of  transferring  energy  from  the  midrange  frequencies  through  wave- 
wave  interactions. 

13.  The  hatched  line  in  Figure  2  represents  the  relocation  of  wave  en¬ 
ergy  between  frequencies  via  the  nonlinear  wave-wave  interactions.  Energy  is 
shifted  in  both  directions  from  region  II.  Energy  shifted  to  higher  frequen¬ 
cies  (region  III)  is  dissipated  through  viscous  and  turbulent  processes.  The 
model  "loses"  this  energy  by  limiting  the  rear  face  of  the  spectrum  to  the 
depth- '.  im i ted  spectral  form.  Energy  shifted  to  the  lower  frequencies  contrib¬ 
utes  to  wave  growth  ltd  the  eventual  shifting  of  the  spectral  peak  to  lower 
frequencies  which  in  turn  allows  the  midrange  to  accept  more  energy  through 
wind  stress.  The  wave-wave  interaction  term  is  sensitive  to  both  t  ,e  energy 
contained  in  the  spectrum  and  the  shape  of  the  spectrum. 

The  model  contains  a  procedure  for  the  transformation  of  swell  wave 
energy.  This  situation  can  arise  either  by  the  input  of  a  bimodal  spectrum  at 
the  onset  of  a  simulation,  by  turning  off  the  winds  at  some  point  during  the 
simulation,  or  by  migration  of  the  spectral  peak  out  of  the  growth  range.  Af¬ 
ter  eacn  Lime-step  the  newly  estimated  wind  sea  spectrum  at  each  grid  point  is 
subtracted  from  the  total  spectrum.  Any  remaining  energy  at  frequencies  less 
than  th"  peak  of  the  wind  sea  spectrum  is  treated  as  swell.  The  model  then 
de" -c-m  i res  the  oeuK  frequency  of  the  swei,  and  the  totai  swell  energy.  It 
compute.-  the  alpha  and  the  wave-wave  interaction  term  for  the  swell,  and  it 
-  b  .  f  •  1--  the  new  p-ak  of  the  swell. 

'A .  ....•  mentioned,  bottom  friction  and  percolation  are  not  considered 

import mt  ; n  the  wind  sea  components  in  this  model;  however,  source  terms  for 
pot  torn  :  riot  ion  and  percolation  are  present,  which  allow  the  user  to  include 
ihi.,*-;  effects  tor  Swell  and  to  select  values  for  the  coefficient  of  friction 
ar'd  ‘.he  coefficient  of  percolation.  The  effect  of  bottom  friction  and  perco- 
i  it  ion  is  to  decrease  "he  energy  levels.  Bottom  friction  and  percolation  are 
no*  coi. -aired  m  iho  wind  sea  oa lance  because  the  empirical  TMA  relationships 
for  wave  grow*  h  ».-e  bei  iev*-d  to  have  accounted  for  their  effect. 

Bound  i pjr_  condi  t  i  or- 

It.  The  ongiria.  version  of  SHALWV  required  that  all  wave  energy  oe 


"grown"  from  an  initially  flat  sea  surface.  This  made  modeling  of  unlimited 
fetch  regions  very  difficult  and  costly,  and  simulations  of  bimodal  spectra 
and  many  prototype  events  were  not  possible.  This  drawback  has  been  rectified 
by  a  feature  in  the  model  which  allows  it  to  initialize  each  grid  point  with  a 
specified  directional  wave  energy  spectrum.  Additionally,  the  incoming  wave 
energy  at  the  boundaries  can  be  updated  during  model  execution  to  simulate  ac¬ 
tual  prototype  boundary  conditions  and  to  prevent  "leakage"  of  wave  energy  out 
the  side  boundaries  (in  an  open  coast  situation).  A  separate  program  which 
provides  the  optional  initial  and  boundary  conditions  for  SHALWV  has  been 
coded  by  Dr.  Robert  E.  Jensen  of  CERC. 

Model  Applications  and  Limitations 


17.  SHALWV  in  its  present  form  is  one  of  the  most  versatile  of  CERC's 
shallow-water  wave  models.  It  is  particularly  well  suited  for  situations  de¬ 
scribed  by  combinations  of  the  following  features  in  finite  depth  water: 

a.  Enclosed  bodies  of  water. 

b.  Open  coasts  (with  additional  codes). 

c.  Arbitrary  bathymetry. 

d.  Nonsteady-state  wind  or  wave  conditions. 

e.  Shallow-water  or  deepwater  wave  growth. 

f.  Swell  wave  transformation. 

g.  Fetch-  or  duration-limited  wave  growth. 

h.  Unlimited  fetch  and/or  unlimited  duration  cases. 

L  Very  lengthy  time  simulations  (hindcasts). 

J.  Rapidly  turning  winds  (e.g.  hurricanes). 

18.  Limitations  of  the  SHALWV  numerical  model  can  be  given  in  two  cate¬ 
gories.  Physical  limitations  refer  to  those  physical  processes  or  situations 
which  are  not  represented  in  the  model.  Some  physical  limitations  are  the 

fol lowing: 

a.  No  interactions  with  currents  (e.g.  vicinity  of  inlets). 

b.  No  interaction  with  coastal  structures. 

c.  No  diffraction  (e.g.  diffraction  behind  islands,  though  islands 
can  be  handled). 

d.  No  variation  of  still  water  level  with  time  (e.g.  storm  surge) 
which  would  affect  the  results  in  very  shallow  water. 


Some  of  the  above  limitations  may  be  rectified  by  future  code  enhancements. 
Computational  limitations  are  functions  of  the  numeric  scheme  employed  in  the 
model,  or  they  can  be  limitations  resulting  from  large  core  memory  require¬ 
ments.  Some  computational  limitations  of  the  present  version  of  SHALWV  are 
that  it 

a.  Allows  no  variable  grids  or  boundary  fitted  grids. 

b.  Requires  a  large,  fast  computer. 

c.  Runs  as  a  batch  Job. 

d.  Has  stability  dependent  upon  Courant  number  (i.e.,  smaller  grid 
means  smaller  time-step  which  means  more  expense). 
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PART  II:  MODEL  EXECUTION 


Program  Organization 


Run  types 

19.  An  important  feature  of  the  present  model  formulation  is  that  all 
time  invariant  interpolation  constants  are  computed  and  stored  in  arrays  prior 
to  the  initiation  of  the  time  loop.  This  is  termed  "preprocessing,"  and  it  is 
an  efficient  scheme  for  time-dependent  models,  but  it  also  rules  out  any  simu¬ 
lations  involving  a  time  varying  water  level.  A  change  of  water  level  would 
require  the  recalculation  of  all  interpolation  constants. 

20.  There  are  three  versions  of  the  JCL  batch  Jobs  which  can  be  run. 

The  first,  PREPJCL,  is  a  preprocessing  run  only.  It  gathers  the  necessary  in¬ 

put  files  and  calculates  only  the  time  invariant  constants  for  the  given  ba¬ 
thymetry  and  specified  parameters.  The  calculated  information  is  saved  as  a 
binary  file  for  use  in  a  production  run.  The  second  batch  Job,  PRODJCL,  reads 
the  file  previously  saved  by  PREPJCL  and  then  proceeds  with  the  wave  simula¬ 
tion.  None  of  the  parameters  set  during  the  preprocessing  run  can  be  reset 
during  a  production  run.  Two  advantages  of  this  arrangement  are  (a)  several 
different  simulations  at  the  same  site  can  be  performed  without  having  to  re¬ 
calculate  the  time  invariant  portions,  and  (b)  the  initial  debugging  of  the 
input  files  is  less  expensive.  The  third  batch  Job,  PRPOJCL,  preprocesses  and 

then  continues  with  the  simulation.  This  Job  stream  is  useful  when  small 

changes  are  made  to  the  input  which  necessitate  preprocessing.  An  example 
would  be  a  change  in  the  bathymetry  file.  Note,  however,  that  PRPOJCL  does 
not  save  the  preprocessing  file  for  future  use.  These  three  run  options  of 
SHAl.WV  are  shown  on  the  simple  flowchart  in  Figure  3.  When  setting  up  runs, 
care  must  be  taken  with  regard  to  some  statements  in  the  JCL  files  that  purge 
old  files,  override  existing  files,  or  initialize  new  files.  Also,  long  simu¬ 
lations  over  the  same  grid  can  create  some  very  large  output  files  that  build 
up  storage  costs. 

Preprocessing 

21.  The  first  step  of  modeL  preprocessing  is  to  read  the  model  input 
parameters.  This  includes  parameters  such  as  numbers  of  rows  and  columns  in 
the  grid,  the  grid  dimensions,  the  time-step,  etc.  There  are  quite  a  few 
parameters  which  allow  the  user  to  select  various  program  options.  At  the 
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Figure  3. 

Model  flowchart 

same  time  the  bathymetry  is  read  and  converted  to  meters.  Reading  of  the  in¬ 
put  data  is  foLlowed  by  the  second  step,  initialization  of  the  working  arrays. 
The  third  step  is  the  calculation  of  the  time-invariant  data  used  by  the  prop¬ 
agation  routines  which  calculate  such  depth  related  parameters  as  phase  and 
group  velocities  for  all  frequencies  at  all  water  grid  points.  They  also  cal¬ 
culate  a  host  of  interpolation  constants  related  to  refraction  and  shoaling. 
These  steps  are  illustrated  in  the  flowchart  in  Figure  3. 

Wave  simulation 

22.  After  preprocessing,  the  program  can  read  optionally  an  initial 
condition  file.  This  file  will  preset  directional  spectral  energy  levels  at 
every  water  grid  point  in  the  model.  This  ability  is  one  of  the  most  useful 
features  of  the  model  because  it  allows  simulations  over  regions  where  the 
boundary  or  initial  conditions  are  known.  An  example  is  an  open  body  of 
water.  Without  this  feature,  the  model  grid  would  have  to  be  big  enough  to 
extend  past  fetch-limited  growth  and  thus  could  become  very  expensive. 

23.  The  program  now  enters  into  the  time-dependent  nested  loops.  In 
the  outer  loop,  the  first  and  second  wind  files  are  read.  This  enables  linear 
interpolation  for  wind  speed  and  wind  direction  at  each  time-step.  Next,  a 
boundary  input  is  read  if  this  option  is  in  effect.  The  program  then  enters 
the  inner  loop  where  the  actual  spectral  calculations  are  performed  for  each 
time-step.  A  series  of  routines  in  this  inner  loop  propagate,  refract,  and 
shoal  wave  energy  and  update  the  energy  as  determined  from  the  source  term 
calculations.  At  the  end  of  each  time-step,  the  one-  and  two-dimensional 
spectra  and  the  one-line  summaries  are  written  to  files  (if  the  user  selects 
this  option);  and  a  new  time-step  is  begun.  When  it  is  time  to  read  another 
wind  file,  the  program  goes  to  the  outer  loop.  Another  enhancement  allows 
output  of  directional  spectra  at  specified  grid  points  to  serve  as  boundary 
input  for  a  nested  grid  application;  however,  running  of  a  nested  grid  re¬ 
quires  an  additional  code.  Finally,  after  all  wind  files  have  been  read  and 
processed,  the  program  will  save  (if  requested)  all  necessary  data  for  a  fu¬ 
ture  warm  start. 

Grid  conventions 

24.  The  grid  used  by  SHALWV  is  a  framework  of  orthogonal  lines  forming 
squares  of  uniform  size.  Any  point  in  the  grid  may  be  denoted  by  I , J  coordi¬ 
nates  as  shown  in  Figure  4.  A  potential  point  for  confusion  in  this  grid  con¬ 
vention  is  the  numbering  of  the  rows  from  the  bottom  to  the  top  of  the  grid. 
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Figure  4.  Model  grid  and  angle  conventions 


Figure  4  also  shows  the  model's  convention  for  directions.  Input  w*nd  direc¬ 
tions  and  output  wave  directions  are  in  a  "going  to"  form  (i.e.,  wind  and  wave 
directions  are  the  direction  from  which  they  came)  with  0  deg  being  parallel 
to  the  rows  with  the  angle  increasing  in  a  counterclockwise  rotation.  For 
example,  in  the  grid  shown,  winds  blowing  from  the  bottom  of  the  grid  to  the 
top  of  the  grid  would  have  an  angle  of  90  deg.  This  direction  convention  is 
also  a  source  of  possible  confusion  since  it  is  contrary  to  the  conventional 
compass  rose.  Consequently,  caution  should  be  exercised  when  (a)  specifying 
grid  points  by  I , J  coordinates,  (b)  specifying  input  wind  angles,  and  (c)  in¬ 
terpreting  output  wave  angles. 

The  parameter  statement 

?5.  SHALWV  contains  a  parameter  statement  in  the  source  listing  which 
acts  to  set  the  dimensions  for  most  of  the  model's  arrays.  An  example  state¬ 
ment  is  given  below. 

PARAMETER  ( IDMN= 16, JDMN= 16, IF =20, 1 A =20, NBPS= 35) 

Each  item  in  the  PARAMETER  statement  is  explained  below. 


Parameter  Name 

Explanation 

I DMN 

Maximum 

number  of  grid  columns 

JDMN 

Maximum 

number  of  grid  rows 

IF 

Maximum 

number  of  frequency  bands 

IA 

Maximum 

number  of  angle  bands 

NBPS 

Max imum 

number  of  input  boundary  points 

14 


The  above  values  must  be  greater  than  or  equal  to  the  corresponding  parameters 
specified  in  the  input  file.  If  a  grid  larger  than  allowed  by  the  parameter 
statement  is  desired,  it  is  necessary  to  change  the  value  in  the  parameter 
statement  to  accommodate  the  increased  size,  and  usually  it  will  be  necessary 
to  make  changes  to  some  of  the  JCL  statements.  The  PARAMETER  statement  occurs 
25  times  in  the  source  code,  and  all  25  statements  must  be  changed  so  that 
they  agree.  When  the  user  fails  to  make  the  values  in  the  PARAMETER  state¬ 
ments  large  enough  for  the  grid,  the  usual  fatal  error  message  is  "indefinite 
result."  This  stems  from  attempting  an  algebraic  operation  with  an  array  ele¬ 
ment  that  does  not  exist.  If  the  values  in  the  PARAMETER  statements  get  too 
large  (e.g.,  a  100  by  100  grid),  even  the  supercomputers  will  not  have  enough 
core  memory  to  run  the  Job. 


Model  Input  Files 


26.  SHALWV  accepts  up  to  a  total  of  six  input  files;  three  originate 
from  the  front-end  computer,  and  three  are  created  and  stored  on  the  Cyber  205 
(or  similar  class  large  computer).  The  input  files  are  given  below. 

Tape  20 

27.  This  file  contains  the  parameters  and  options  selected  by  the  user. 
It  is  required  for  all  three  batch  Job  streams.  It  is  a  front-end  file. 

Tape  25 

28.  This  file  contains  the  bathymetry.  It  is  required  for  batch  jobs 
PREPJCL  and  PRP0JCL.  It  is  a  front-end  file. 

Tape  21 

29.  This  file  specifies  the  winds  over  the  grid.  It  is  required  for 
batch  Jobs  PR0DJCI.  and  PRP0JCL.  It  is  a  front-end  file. 

Tape  33  (Prep) 

30.  This  is  a  binary  file  created  by  the  batch  job  PREPJCL  and  stored 
on  the  large  computer.  It  contains  all  the  preprocessing  information,  and  it 
is  called  by  the  batch  Job  PR0DJCL.  This  file  will  not  be  created  by  PRP0JCL. 
Tape  15  (Warm) 

31.  This  is  a  binary  tape  which  can  be  written  by  either  PRODJCL  or 
PRP0JCL.  It  contains  information  necessary  for  a  "warm  start,"  i.e.,  the  con¬ 
tinuation  of  a  simulation.  This  is  saved  on  the  large  computer. 
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Tape  31  (  Bndy  in  ) 

32.  This  is  a  file  containing  specified  initial  and  boundary  condi¬ 
tions.  It  is  saved  on  the  large  computer,  and  it  can  be  called  by  either 
PRODJCL  or  PKPOJCL.  This  file  is  created  by  a  separate  program. 

Model  parameters  input  (Tape  20) 


33.  Input  for  preprocessing,  PREPJCL.  The  following  input  parameters 
need  to  be  specified  for  a  preprocessing  run.  Figure  5  is  a  sample  input  file 
used  for  model  runs  at  the  Field  Research  Facility  (FRF)  at  Duck,  North  Caro¬ 
lina.  Each  parameter  is  explained  below.  Values  given  in  parentheses  are 
from  the  example  in  Figure  5  and  are  included  only  for  clarification.  They 
are  not  default  values. 


14  15  16  14  0  1 

3.  10.  200.  0.3048 

0.060  0.010 
00000000000000 
01111111111110 
nu  in  mi  u  io 
0 1 1 1 1 1 1 1 1 1 1 1 10 
01111111111110 
01111111111110 
oi  n  mi  ii  u  in 
omiiiinuio 
oiiim  iiiuio 
0 1 1 1  ill 1 1 1 1  n  o 
01111111111110 
01111111111  110 
01  111  1 1  111  1  1  10 
01 1 1 1 1 1 1111 1 10 
00000000000000 

Figure  5.  Sample  parameter  input 
file  for  preprocessing  run 


J_st_  Car'd  -  Format!  14  )  . 

1  PREP  -  Run  type  (1),  IPPF.P  =  1  preprocess  and  stop 

(this  example) 

=  2  production  run 
=  3  preprocess  and  run 

c’nq  ml  -  Format  ( 6  1 4  )  . 

NX  -  Number  of  columns  in  grid  (14) 

NY  -  Number  of  rows  in  grid  (15) 

NANG  -  Number  of  directional  angle  bands  (16) 
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NFRC  -  Number  of  frequency  bands  (14) 

INDELF  -  Flag  on  frequency  band  read  (0) 

=  0  internally  computes  frequency  bands 
=  1  reads  frequency  bands 

IDEPD  -  Flag  on  water  depth  read  of  Tape  25  ( 1 ) 

=  0  read  in  constant  water  depth 
=  1  read  in  variable  water  depths 

c.  3rd  Card  -  Format( 3F8. 0,F8 .4 ) . 

DL  -  Distance  between  grid  points  in  kilometers  (3-0) 

XN  -  A  scale  factor  used  in  refraction  computations  ex¬ 
pressed  as  a  percentage  of  DL  (10.0)  (A  value  of  10  is 
recommended.  A  larger  number  increases  computational 
costs;  a  lower  number  decreases  refraction  computation 
accuracy . ) 

DT  -  Time-step  in  seconds  (200.0) 

DEPFAC  -  Factor  to  convert  input  depths  to  meters  (0.3048) 

d.  4th  Card. 

IF  INDELF  ::  0,  Format(2F6.3) 

0MST  -  Starting  frequency  band  (0.060) 

DOMEG  -  Frequency  interval  (0.010) 

IF  INDELF  =  1,  Format( 10F6. 3) 

0MEG(K),  K=1,NFRC  -  Read  in  frequency  bands,  can  take  two  cards 

e.  5th  Card  -  Format(21I1). 

I B0UND( I , J )  -  Array  designating  grid  points  as  water,  land,  or 
boundary  points  (All  outer  grid  points  must  be 
denoted  as  boundary  points.  This  array  must  be 
the  same  size  as  the  bathymetry  array.  ) 

=  0  boundary  point 
c  1  water  point 

This  ends  the  Tape  20  input  file  for  the  batch  Job  PREPJCL. 

34.  Input  for  Production  Run,  PRODJCL.  The  following  input  parameters 
need  to  be  specified  for  a  production  run.  Figure  6  is  a  sample  input  file 
used  for  model  runs  at  the  FRF.  Each  parameter  is  explained  below.  Values 
given  in  parentheses  are  from  the  example  in  Figure  6  and  are  included  only 
for  clarification.  They  are  not  default  values, 
a .  1st  Card  -  Format! 14). 

IPREP  -  Run  type  (2),  TPREP  =  1  preprocess  and  stop 

=  2  production  run  (this  example) 

=  3  preprocess  and  run 
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r  2  cold  start,  save 

-  3  warm  start,  save 

IBTY  -  Flag  for  boundary  condition  input  and  output  (2) 
=  0  no  input  or  output  boundary  conditions 
=  1  output  boundary  information  (Tape  30) 

=  2  input  boundary  information  (Tape  31) 

=  3  input  and  output  boundary  information 

NORD  -  Flag  for  wind  input  read  (2) 

=  1  constant  over  time  and  space 
=  2  variable  over  time,  constant  over  space 
=  3  variable  over  time  and  space 

NOBE  -  Flag  for  bottom  effects  (1) 

=  1  no  bottom  effects 

-  2  include  bottom  effects 

3rd  Card  (optionaL)  -  Format(2F10.5) . 


IF  NOBE  -  2,  this  card  is  read;  otherwise  do  not  include  this 
card . 

BCOF  -  Bottom  friction  coefficient 

PCOF  -  Bottom  percolation  factor  in  cm/sec 

4th  Card  -  Format(I4). 

NTMS  -  Number  of  time-steps  between  output  results  (18) 

NH -  Number  of  hours  between  input  winds  (1) 

MXHR  -  Maximum  number  of  wind  inputs  (10) 

MSTA  -  Number  of  special  output  locations,  not  to  exceed 
10  (3)  (These  are  the  grid  points  where  detailed 
spectral  wave  data  are  saved.) 

5th  Card(s)  -  Format(2I4). 

1 OUT ( K )  -  Column  (I)  location  of  special  output  point 

JOUT(K)  -  Row  (J)  location  of  special  output  point 

Them  must  be  a  total  of  MSTA  of  these  cards,  one  for  each 
location . 

Note:  The  cards  from  this  point  pertain  to  the  input  or  output 
of  boundary  information.  To  make  use  of  the  data  specified  by 
these  cards  additional  codes  are  reouired. 

6th  Card  (optional)  -  Format(I4). 

IF  IBTY  :  1  or  3,  this  card  is  read;  otherwise  do  not  include 
this  card  (not  included  in  Figure  6  example). 

NOBS  -  Number  of  output  locations 

7th  Card(s)  (optional)  -  Format(2I4). 

IF  IBTY  =  1  or  3,  a  total  of  NOBS  cards  must  be  provided,  one 
card  per  output  boundary  point. 

L0UTB(K,1)  -  Column  (I)  location  of  boundary  output  location 
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L0UTB(K,2)  -  Row  (J)  location  of  boundary  output  location 
where  K  =  1,  NOBS 

h.  8th  Card  (optional)  -  Format(l4). 

IF  IBTY  :  2  or  3,  this  card  is  read;  otherwise  do  not  include 
this  card  (included  in  this  example). 

NIBS  -  Number  of  input  boundary  locations  (35) 

9th  Card(s)  (optional)  -  Format(2I4). 

IF  IBTY  :  2  or  3,  a  total  of  NIBS  cards  must  be  provided,  one 
card  per  each  input  boundary  point. 

LINB(K,1)  -  Column  (I)  location  of  boundary  input  location 

LINB(K,2)  -  Row  (J)  location  of  boundary  input  location 

where  K  =  1 ,  NIBS 

CAUTION:  these  cards  must  be  specified  in  the  exact  order 
specified  in  the  boundary  generator  program. 

This  ends  the  Tape  20  input  file  for  the  batch  job  PRODJCL. 

35.  Input  for  preprocessing  and  production  run,  PRPOJCL,  The  input 
file  for  a  preprocessing  and  production  run  is  a  combination  of  the  prepro¬ 
cessing  cards  and  the  production  cards.  No  example  is  provided  for  this  case, 

a .  1st  Card  -  Format(l4). 

I  PREP  -  run  type  (2),  I  PREP  =  1  preprocess  and  stop 

=  2  production  run 
=  3  preprocess  and  run  (this 
example) 


b.  Cards  2  through  5. 

These  cards  will  be  identical  to  cards  2  through  5  of  a  pre¬ 
processing  ( PREPJCL )  run. 

c .  Cards  6  through  13. 

These  cards  will  be  identical  to  cards  2  through  9  of  a  produc¬ 
tion  (PRODJCL)  run. 
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This  ends  the  Tape  20  input  file  for  the  batch  job  PRPOJCL. 

Bathymetry  input  (Tape  25) 

36.  The  bathymetry  grid  is  an  array  containing  positive  values  of  depth 
at  each  grid  point.  Depths  can  be  in  any  units,  but  units  cannot  be  mixed 
within  the  grid.  Large  grids  are  handled  by  reading  successive  cards  until 
the  entire  row  is  read.  The  next  row  begins  on  a  new  card.  An  example  is 
3hown  in  Figure  7,  which  is  the  bathymetry  file  (in  feet)  used  for  simulations 
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0.0 

35.0 
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55.0 

65.0 

60.0 

84.0 

81. 0 
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56.0 

72.0 
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78.0 

80.0 

78.0 

67.0 

65.0 

76.0 

80.0 

82.0 

84.0 

65.0 

71.0 

77.0 

85.0 

84.0 

87.0 

70.0 

79.0 

85.0 

83.0 

85.0 

90.0 

60.0 

70.0 

92.0 

83.0 

85.0 

92.0 

62.0 

80.0 

87.0 

84.0 

85.0 

90.0 

65.0 

60.0 

84.0 

82.0 

84.0 

92.0 

61.0 

79.0 

83.0 

80.0 

85.0 

84.0 

65.0 

64.0 

80.0 

77.0 

80.0 

87.0 

68.0 

74.0 

71.0 

78.0 

78.0 

78.0 

62.0 

78.0 

80.0 

85.0 

90.0 

79.0 

61.0 

63.0 

82.0 

70.0 

90.0 

70.0 

61.0 

63.0 

82.0 

70.0 

90.0 

70.0 

Figure  7.  Sample  bathymetry  grid 
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File  input 


IF  IDEPD  = 
DEPCO 
IF  IDEPD  - 
DEP(I,J)  - 


0  (Tape  20),  Format(F6.1) 

Uniform  (constant)  water  depth  over  entire  grid 
1  (Tape  20),  Format ( 20F6 . 1 ) 

Array  containing  water  depths  at  all  grid  points  (This 
array  must  be  the  same  size  as  the  IBOUND  array  on 
Tape  20. ) 


Wind  input  (Tape  21) 

37.  The  wind  input  file  is  required  for  all  runs  of  the  batch  jobs 
PHODJCL  and  PRPOJCL..  A  wind  field  can  be  input  as  either  a  single  card  spec¬ 
ifying  a  uniform  and  constant  wind  speed  and  wind  direction  over  the  entire 
grid,  or  it  can  be  in  the  form  of  arrays  which  give  the  wind  speed  and  direc¬ 
tion  at  each  grid  point.  The  wind  file  can  contain  any  number  of  wind  fields, 
but  it  must  contain  at  Least  two.  When  wind  fields  are  not  input  at  every 


time-stop,  a  linear  interpolation  is  performed  to  obtain  the  intermediate  val¬ 
ues  at  each  grid  point. 


a .  Uniform  wind  fields  -  Format ( 2F6 ,  1 ) . 

lor  values  of  NORD  :  1  or  2  (from  Tape  20),  wind  speed  and  wind 

direction  are  assumed  not  to  vary  over  the 
grid.  The  input  15  then 

TSWS  -  Urnfcir.  wind  sro»*.i  in  km  •  j  cv--r  u-nt : :v  f*ric  assumed 
to  b*;  ovorwater.  neutrally  stable  atrnsphoru ?  cjr.u_- 
t.  ior.s  at  the  ‘  U-rr.  •  I  ova  tier. 

TSDIR  -  Uniform  wind  direction  m  degrees  over  entire  grid  ex¬ 
pressed  relative  to  direction  convention  discussed 
before  J 


IF  NORD  -  1,  only  one  card  is  read. 


m 


"  v v-w vr»  v*  * ."svr,Tv^€ni  \p«  VV\>  ^  .  J  ^  .  ■-' 


IF  WORD  =  2,  there  will  need  to  be  MXHR  cards  with  the 
above  format. 

Figure  8  is  an  example  of  the  wind  file  for  uniform 
wind  fields  changing  every  hour  for  9  hr. 

b.  Variable  wind  fields.  When  WORD  =  3  (from  Tape  20), 
the  wind  fields  are  assumed  to  be  variable  over  space 
and  time.  Each  wind  field  must  be  organized  as  follows: 

1st  Card  -  Format(IIO) 

ID  -  date/time  code 

Wind  Speed  Array  -  Format(2 1F6 . 1  ) 

WS(I,J)  -  wind  speeds  in  knots  over  the  entire  grid 

Wind  Direction  Array  -  Fermat ( 21 Fo.1 ) 

WD(I,J)  -  wind  direction  in  degrees  over  the  entire 
grid  expressed  relative  to  direction  convention  dis¬ 
cussed  before. 

Figure  9  gives  an  example  of  a  variable  wind  field.  The  wind 
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grids  with  more  than  21  columns,  the  first  21  columns  of  ail 
rows  are  read  first.  This  is  followed  by  reading  columns  22-42 
of  all  rows,  etc.  The  arrays  in  the  wind  field  must  agree  in 
dimension  with  the  bathymetry  array.  The  J=1  row  is  the  bottom 
row,  just  as  in  the  bathymetry  grid.  There  must  be  MXHF)  (Tape 
20  parameter)  wind  fields  present  in  the  wind  file. 


Model  Output 


38.  SHALWV  outputs  up  to  six  files;  three  are  always  output  to  the 
front  end  computer  and  saved  under  the  user's  directory,  one  is  called  back 
interactively  from  the  Cyber  205,  and  two  optional  output  files  reside  on  the 
205  system  for  future  use.  The  output  files  are  as  follows: 


File  Name 
Output 
Tape  10 
Tape  1 1 
Tape  13 
Tape  30 

Tape  15 


_ Information  It  Contains _ 

Computer  run  and  the  input  data 
Two-dimensional  energy  spectra 
Statistical  summaries 
One-dimensional  energy  spectra 

Two-dimensional  spectra  meant  to  serve  as  boundary 
input  for  a  future  run  (optional) 

Necessary  information  for  a  warm  start  (optional) 


Mod el  output  til e 

39.  This  file  must  be  retrieved  interactively  after  the  run  has  been 
completed.  Its  primary  contents  are  (a)  the  compiled  (unmapped)  source 
listing,  (b)  verification  of  program  options  entered  by  the  parameter  input 
file,  (c)  various  statements  indicating  the  progress  of  the  program  during 
execution,  (d)  tables  of  the  combined  significant  wave  height,  peak  spectral 
wave  period,  and  mean  wave  direction  at  every  grid  point  every  time  a  new  wind 
file  is  input  during  execution,  and  (e)  a  dayfile  giving  the  sequence  of  job 
execution  on  the  Cyber  205  and  on  the  front-end  machine.  Only  items  a,  b,  and 
e  are  included  in  a  preprocessing  run. 

40.  Figure  10  shows  an  example  of  the  start  of  the  output  directly  below 
the  compiled  listing.  Figure  11  is  a  sample  of  the  periodic  output  over  the 
entire  grid. 

41.  After  this  file  has  been  retrieved  interactively,  it  must  be  saved 
in  the  user's  directory  if  so  desired.  It  is  not  automatically  placed  in  the 
directory.  Significant  storage  savings  can  be  made  by  deleting  the  compiled 
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THIS  RUN  CONSISTS  OF  A  READ  PREPROCESS  AND  RUN  WAVES  . 
**  NOTE  ALL  INFORMATION  READ  FROM  FILE. 

THIS  RUN  HAS  A  COLD  START  WITH  A  SAVE. 

NUMBER  OF  ROWS  IN  GRID  -  15 

NUMBER  OF  COLUMNS  IN  GRID  -  14 

DISTANCE  IN  METERS  BETWEEN  GRID  POINTS  -  3000.00000 

THERE  ARE  14  FREQUENCIES 

FREQUENCY  BANDS  WERE  COMPUTED  INTERNALLY. 


.0600  .0700  .0800  .0900 

.1000 

.1100 

.1200 

.1300 

.1400 

.1500 

.1600  .1700  .1800  .1900 

THE  14  FREQUENCY  BAND  WIDTHS 
.0100  .0100  .0100  .0100 

ARE 

.0100 

.0100 

.0100 

.0100 

.0100 

.0100 

.0100  .0100  .0100  .0100 
THERE  ARE  16  DIRECTION  BANDS. 

EACH  DIRECTION  BAND  WIDTH  IS  22.5  DEGREES. 

.0  22.5  45.0  67.5  90.0  112.5  135.0  157.5  180.0  202.5 

225.0  247.5  270.0  292.5  315.0  337.5 
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Figure  10.  Beginning  of  a  sample  output  file 
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Figure  11.  Sample  of  periodic 
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source  code  listing  from  the  output  file. 


Two-d imens iona 1  spectra  (Tape  10) 

4?.  This  file  contains  energy  density  values  (units  of  m-m-sec/rad)  for 
two-dimensional  spectra  at  each  special  output  location  for  each  specified  time 
increment.  In  the  example  shown  in  Figure  12,  the  first  line  is  the  date/time 
code  and  the  I,J  location  of  the  special  output  point.  This  example  two- 
dimensional  spectrum  has  16  angle  bands  (columns)  beginning  with  0  deg  and 
incremented  by  22.5  deg.  The  14  frequency  bands  (rows)  begin  with  0.1  Hz  and 
are  incremented  by  0.01  Hz.  This  file  is  automatically  saved  as  an  indirect 
access  file.  This  file  is  best  displayed  or  printed  using  132  column 
terminals. 
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Figure  12.  Sample  two-d  i 
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St  it  i  st  i cal  Summaries  ( Tape  1 1 ) 

4j .  This  file  contains  one-line  summaries  of  significant  wave  parameters 
(e.g.  significant  wave  height,  HMO;  peak  spectral  period,  TP:  and  mean  direc¬ 
tion,  TH)  for  both  sea  and  swell  at  each  special  output  location  for  each 
specified  time  increment.  An  example  from  Duck,  North  Carolina,  is  given  in 
Figure  13.  All  one-line  summaries  for  a  given  time-step  are  listed  by  output 
location  I , J .  These  sets  of  lines  are  repeated  in  chronological  order.  Each 
line  also  gives  the  wind  speed,  WS,  ana  the  wind  direction,  V.D,  at  that  grid 
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13.  Sample 

summary 

file 

output 

point  at  that  time-step.  This  file  is  automatically  saved  as  an  indirect  ac¬ 
cess  file. 

One-dimensional  spectra  (Tape  13)  ' 

44.  This  file  contains  energy  density  values  (units  of  m-m-sec)  for  one 
dimensional  spectra  at  each  special  output  location  for  each  specified  time 
increment.  In  the  example  given  in  Figure  14,  the  first  line  gives  the  water 
depth  in  meters,  the  wind  speed  in  knots,  and  the  wind  direction  (in  5-deg 
increment  bands  with  band  number  1  centered  on  0  deg).  The  second  line  is  the 
date/time  code  and  the  I , J  location  of  the  special  output  point.  The  next  two 
lines  give  the  energy  density  of  the  14  frequency  bands  beginning  at  0.1  Hz 
and  incremented  by  0.01  Hz.  This  file  is  automatically  saved  as  an  indirect 
access  file. 

Boundary  output  file  (Tape  30) 

45.  This  optional  file  contains  two-dimensional  spectra  at  specified 
time-steps.  This  output  file  (after  manipulation  by  another  code)  Decomes  the 
boundary  input  file  (Tape  31)  for  a  smaller  nested  grid  run  which  also  re¬ 
quires  another  code. 

Warm  start  file  (Tape  15) 


46.  This  output  file  is  described  in  the  section  on  model  input  files. 
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Figure  14.  Sample  one-dimensional  spectra  output 
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APPENDIX  A:  SHALWV  CYBERNET  EXECUTION 


1.  This  Appendix  gives  the  details  of  setting  up  and  running  the  SHALWV 
model  on  the  Cybernet  Computer  System.  Most  of  this  material  is  specific  to 
the  Cybernet  system. 


Job  Control  Language  Files 

2.  The  Job  control  language  (JCL)  file  is  a  batch-run  file  that  gathers 
up  the  necessary  input  files,  sends  them  off  to  the  Cyber  205  for  execution, 
and  finally  retrieves  and  saves  the  generated  output  files.  The  preprocessing 
JCL  ( PREPJCL)  is  shown  in  Figure  A1,  production  JCL  (PRODJCL)  is  listed  in 
Figure  A2,  and  the  preprocess  and  run  JCL  (PRPOJCL)  is  given  in  Figure  A3. 

3.  In  the  JCL  files  listed  in  Figures  A1,  A2,  and  A3,  the  sequence 
AF=CERCFIL  appears  several  times.  This  is  a  reserved  file  namt  in  the  link- 
pool  that  the  user  must  create  prior  to  running  the  model. 

4.  Further  explanation  of  the  JCL  statements  can  be  found  in  the  Cyber¬ 
net  system  documentation.  Be  aware  that  JCL  commands  and  conventions  are  oc¬ 
casionally  changed,  so  the  example  JCL  files  may  be  outdated. 


PREPJCL. 

USER(  U-xxxxxx  , PA-yyyyyy)ADY. 

RESOURCE , JCAT-P4 , TL-900 . 

CHARGE,  CKROF.CC.xyzxyz. 

PURGE, PREP. 

U  NK ,GET( CERC-MODEL/ FM-KOK , AF-CERCFI L , DD-C6 ) 
LINK,GKT(TAPE20-PREPIN/FM-K0E,AF-CERCFTL,DD-C6) 
LINK,GET(  TAPE2  5-DUCKL/  FM-KOE ,  AF-CERCF IL ,  DD-C6  ) 
REQUEST , PREP/9999 , RT-W. 

FORTRAN, l-CERC,B-CERCfi/#50. 

LOAD , CERCB , CN-CERCGO/ 5000 .CDF-4000 .GRLPALL-  . 
CERCGO.TAPE33-PRFP. 

DEFINE.PREP/9999, RT-W. 

SUMMARY. 

EXIT. 


Figure  A1.  Preprocessing  JCL  (PREPJCL) 


PROD JCL. 

USKR(  [J-xxxxxx ,  PA-yyyyyy)ADY. 

RESOURCE  ,  .TCAT-P4  ,  TL-900 . 

CHARGE, CEROEGC.xyzxyz. 

LI  NK ,  GF.T( CERC-MODEL/ FM-KOE , AF-CERCF I L , DD-C6 ) 

LI NK,GET< TAPE20-PRODIN/FM-KOE, AF-CERCFIL, DD-C6) 

L  INK, GF.T(TAPF,21-WlNn3/FM-KOE,  AF-CERCFIL,  DD-C6) 
ATTACH, PREP. 

ATTACH, WARM. 

ATTACH, BNDYIN. 

FORTRAN, I-CERC,B-CERCB/#50. 

REQUEST .TAPE  10/ 1000 . 

LOAD ,  CERCB ,  CN-CERCGO/  5000  .GRLPALL-  , CDF-4000 . 
CERCCO.TAPF 3  3-PREP, TAPE  15 -WARM, TAPE31-BNDYIN. 

LINK,  REP  LACE ,  TAPE  1 0-TAPE 1 0/ FM-KOE ,  DD-C6  , AF-CERCF IL 
LINK, REPLACE , TAPE  1 1 -TAPE i i / FM-KOE , DD-C6 , AF-CERCFIL 
LINK, REPLACE, TAPE  1  3-TAPE  1 3/ FM-KOE ,DD«C6 , AF-CERCFIL 
PURGE, WARM. 

SUMMARY. 

EXIT. 

L I  NK ,  REP l.AC K  ,  TAPE  1 0-TAPE 1 0/  FM-KOE ,  DD-C6  ,  AF-CERCF  IL 
LINK, REPLACE, TAPE1 1 -TAPE  1 1 /FM-KOE, DD-C6 , AF-CERCFIL 
LI NK , RE  P LACK , TAPE  1  3-TAPE  1 3/ FM-KOE , DD-C6 , AF-CERCF IL 

Figure  A2.  Production  JCL  (PRODJCL) 


PRPOJCL. 

USER( U-xxxxxx ,PA-yyyyyy) ADY. 

RESOURCE , JCAT-P4 , TL-900 . 

CHARGE , CEROKGO , xyzxyz . 

LINK, GKT(CERC-CKRC2/ FM-KOE, AF-CERCF LL,DD=C6) 

L I  NK,Gr'T(  TAPK20- PR  PRO  IN/ FM-KOE,  AF-CERCFIL,  DD-C6) 
LINK , CET( TAP E 2  S-DUCKL/ FM-KOE , AF-CERCF IL, DD-C6) 

LI SK,GKT(TAPE2 l -W I ND3/ FM-KOE, AF-CERCFIL, DD-C6) 
FORTRAN,  I  -  CERC  ,  B-CERCB/  //50  . 

LOAD , CKRCh , CN-CERCGO/ #5000 .CDF-4000 .GRLPALL-  . 

C.  r,  K  <  ’( 

LINK  ,  Ki.P LACK  , TAPS  1  O-TAPE 10/ FM-KOE  ,  DD— CF  ,  AF-CERCFIL. 
LINK, REPLACE, TAPE! 1 -TAPE l 1 /FM-KOE, DD-C6, AF-CERCFIL. 
LINK,  RE  P  [.AC  K  ,  TAPE  1  3-TAPE  13  /  FM-KOE ,  DD-C6  ,  AF-CERCF  I L . 

SUMMARY. 

EXIT. 

LINK,  RHP  LAC  ,  TAPE  i  0-TAPE  1  0/  FM-KOE  ,  DD-C6  ,  AF-CERCFIL. 
LINE , REPLACE , TAPE  1 l-TAPEl i /FM-KOE , DD-C6 , AF-CERCFIL. 
l.TNK ,  REPLACE  ,  TAPE  1  3-TAPE  1  3  /  FM-KOE  ,  DD-C6  ,  AF-CERCFIL. 


P  igure  .*3.  Preprocess  und  run  JCL  (PRPOJCL) 


i-v  V  V  W*.  A 

i 


Executing  a  JCL  File 


5.  To  interactively  run  this  batch  job,  issue  the  following  commands 
from  the  terminal. 

GET, JCLFILE 
SUBMIT, JCLFILE.T 


Example:  G£T,PREPJCL 

SUBMIT, PREPJCL.T 


To  check  on  the  status  of  the  job  from  an  interactive  terminal,  type: 

LINK, ENQUIRE 

When  the  Job  is  available  for  output  to  the  terminal,  a  seven-character  job 
name  will  be  displayed  by  the  LINK, ENQUIRE  command  (example:  AF5B6B2).  To 
bring  the  output  file  to  the  terminal,  type 

QGET.xyz  where  xyz  are  the  last  three  characters  of  the  file¬ 

name. 

Example:  QGET,6B2 

This  command  retrieves  the  output  file  and  makes  it  local.  Now  it  can  be  made 
an  indirect  access  file  with  the  interactive  command 

REPLACE,  AF5B6B2=myfile  (for  example) 

or  it  can  be  placed  into  the  editor  for  examination 

XEBIT, AF5B6B2  (for  example) 

6.  The  primary  output  data  files  are  automatically  saved  as  indirect 
access  files.  And  each  time  the  Job  is  run,  the  new  output  files  will  over¬ 
write  the  old  files  unless  the  old  files  are  renamed. 


A3 


